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^ '. Abstract 

^O \ We consider an inextensible, semifiexible polymer or worm-like chain which is confined in the 

transverse direction by a parabolic potential and subject to a longitudinal force at the ends, so 
that the polymer is stretched out and backfolding is negligible. Simple analytic expressions for the 
partition function, valid in this regime, are obtained for chains of arbitrary length with a variety of 
S \ boundary conditions at the ends. The spatial distribution of the end points or radial distribution 

'^ ■ function is also analyzed. 
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I. INTRODUCTION 

The statistical properties of biological polymers in channels with diameters comparable 
with the persistence length have been studied in several recent experiments l|-|8| . When the 
channel diameter is much smaller than the persistence length, the polymer is stretched out 
in the channel, so that its longitudinal length is only slightly shorter than its contour length 
and there is negligible backfolding. This is the regime studied theoretically in this paper. 

We consider the simplest model for a confined biopolymer - an inextensible, semiflexible 
filament or worm-like chain with persistence length P and contour length L, confined in 
the transverse direction and subject to a longitudinal stretching force or tension r. If the 
polymer is tightly confined and/or strongly stretched, the line or filament by which we model 
it is almost straight. Each such configuration correspond to a single valued function r(t), 
where (x, y, t) are Cartesian coordinates, and r = (x, y) specifies the transverse displacement 
of the polymer from the t axis. Since the slope v = df/dt satisfies \v\ <^ 1, the difference i^ 
between the contour length L = J^ ^^ dt [1 + v{tY] and the longitudinal length i?|| may be 
replaced by 

i = L-R\\ = ]^l dt v{t)\ (1) 

and the Hamiltonian l-i of the worm-like chain |9i] simplifies to 



kuT 



BJ- Jo 



dt 



2 \dty 2 \dt ^ ' 



(2) 



The three terms on the right side of Eq. ([2]) correspond to the bending energy, the potential 
energy —t{R\\—L) associated with the stretching force r = ksTf, and the confining potential 
per unit length, all divided by ksT. For a polymer in a channel with hard walls, V{r) takes 
the values and oo for r inside and outside the channel, respectively. 

Note that in Eqs. ([T]) and (|2]), we have replaced the upper limit R\\ in the integrals over 
t by the contour length L, as is correct to leading order in the regime \v\ ^ 1, and in this 
approximation identify the integration variable t with the arc length. Taking this same point 
of view, we interpret the path integral Z = J Dx J Dy exp {—'H/kBT) with T-C given by Eq. 
([2]), studied in detail below, as the partition function of a polymer of fixed contour length L 
but fluctuating longitudinal length /2|| = L— ^ f^ dt v{ty, and from Z determine the mean, 
variance, and distribution of R\\. 



The properties of a semi-flexible polymer in cylindrica 
dinal stretching force, have been studied by simulations 
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. channels, but without the longitu- 
lOMlTJI and calculated theoretically 
19| with the channel replaced by a parabolic confining potential. The dynamical evo- 



lution of the confined polymer under sudden changes of a longitudinal stretching force or 
relaxation of the confinement is studied in Ref. 20 1. 



In this paper we assume that the confining potential in the Hamiltonian ([2]) has the 
parabolic form 

V{r) = \ {Kx' + byy') . (3) 

In Section II an exact analytical result is obtained for the corresponding partition function 
Z (r, v; tq, vq] t) for a polymer of contour length t = L, with arbitrary fixed values r, v and 
"rojVo of the position and slope at t and at t = 0, respectively. An analogous result, but 
without the second, longitudinal force term in the Hamiltonian (12]), was derived in Ref. 
18| . Here we consider the complete quadratic Hamiltonian, quadratic in all three quantities 
dr/dt, dr/dt and r. We note that the boundary condition of fixed end points and end 



slopes is more general than the boundary conditions considered in Refs. [17|, |19| , which are 
basically periodic, so that the Hamiltonian ([2]), expressed in terms of Fourier components 
of r*(t), is diagonal. In Section HI compact analytical expressions for the partition function 
for those and several other boundary conditions, including the experimentally relevant case 
of freely fiuctuating ends, are derived from our result for Z (r, v] ro,vo;t). In Section IV we 
study the equilibrium response of the polymer to the longitudinal stretching force and show 
how to calculate the endpoint distribution or radial distribution function by making the 
replacement / — t- / + s in the partition function and then performing the inverse Laplace 
transformation s -^ L — R\\. Section V contains closing remarks. 



II. EVALUATION OF THE PARTITION FUNCTION 



For the Hamiltonian ([2]) with confining potential ([3]), the fiuctuations in the x and y 
directions are statistically independent, and the partition function factors in the form 

Z (f, v; fo, vo; t) = zfl^{x, v^; Xq, v^] t)zfl^{y, Vy] yo, Vyo; t), (4) 

where 



'(2) 



Zj^l {x,v;xo,Vo;t)= / Dx exp 



1 



dt 



pi'^y+f('4X+bx^ 



dt^ 



(5) 



is the corresponding path integral in the two-dimensional space {x,t). Except for the term 



with coefficient / in the action, the path integral ([5]) is the same as considered in Ref. 18 1 
and can be evaluated in the same way. As in the Feynman-Hibbs treatment of the quantum 
harmonic oscillator [2l|, we substitute x{t) = x*{t) + C,{t) in Eq. ([S]), where of all the paths 
with the prescribed end points and end slopes, x*{t) is the "classical" path which minimizes 
the action. Since the action is quadratic in x and its ffist two derivatives and minimized by 
x*{t), and since ^(t) and d^{t)/dt vanish at the end points, 

Z)l{x,v\xo,VQ]t) = Z].;;(0,0;0,0;t)exp[-S'*(x,w;zo,wo,t)], (6) 

where S*{x, v; xq, vq, t) is the action for the classical path. 
The classical path satisfies the Euler-Lagrange equation 



which has the general solution 



where 



\t)=Y,C^g^{t), (8) 



X 

fc=l 



giit) = cosh (at) cos (/3t) , (9) 

^2(t) = cosh (at) sin (/3t) , (10) 

gsit) = sinh (at) cos (/3t) , (11) 

g^it) = sinh (at) sin (/3t) , (12) 

and 

a=y^, /3=/^' l = y{bP)-'/^ i = b^/'p-^/H. (13) 

Choosing the expansion coefficients ci,...,C4 in Eq. (J8]) so that x*(t) reproduces the pre- 
scribed positions and slopes at the ends of the polymer and then evaluating the action, we 



obtain 



where 



and 



S*{x,v;xo,vo,t) = D~^ 

X { [a^^ sinh (2at ) + 13'^ sin (2/3t )] (x^ + xl) 
—2 [a^^sinh^ (at) + /3~^sin^ (/^O] (^^ ~ 2:0^0) 
+ [a-^ sinh {2ai) - f3'^ sin (2/3t )] {v^ + v^) 
—4 [z?""'^ cosh (at ) sin (/3t ) + a~^ sinh (at ) cos (/3t ) ] xxq 
—4 [a^'^(3^^ sinh (at ) sin (/3t )] {xvq — vxq) 
+4 [/3"^ cosh [at) sin (/3t ) - a'^ sinh (at ) cos (/3t )] vvq] , (14) 

L> = 2 [a^2 sinh^ (at) - (5'^ sin^ (/3t )] , (15) 



a; 



53/8pl/8^^ t) = 6l/8p3/^. (16) 



The prefactor Z\l (0, 0; 0, 0; t) in Eq. ([6]) is determined by substituting Eq. ([6]), with S"* 
given by Eqs. ([H]) and flTSj) . in the Fokker-Planck equation 

Solving the resulting differential equation for Zj- ^ (0, 0; 0, 0; t) with the initial condition 

Z)l (x, f ; zo, t^o; 0) = 5{x - Xo)5{v - vq), (18) 

we obtain 

Zfl (0, 0; 0, 0; t) = Ti-\hPf/^ [a"' sinh' (at ) - ^^ sin' (/^O] ~'^' • (19) 

This completes the exact evaluation of the polymer partition function (jlj). The result 
is given in Eqs. (jH]), fll3l) - fll6l) . and fll9p . We have checked with Mathematica that it does 
indeed satisfy the Fokker-Planck equation (TT71) with initial condition (TT8|) . In the limit / — )■ 
it is consistent with Eq. (13) of Ref. [iSJ], and the value of the unspecified normalization 
constant J\f in that equation is found to be A/" = (27r') . 

The function /3 of the variable 7 defined in Eq. ( fT3l) has a square root branch point at 
7 = 4/(6P)^^/' = 1. No singularity in the partition function is expected at this value of 7, 



and that there is none is clear from Eqs. (fT^ . (fT5|) . and (fT9|) . The quantity (3 only occurs 
in the combinations like /3~^ sin (/3t) and cos (/3t). Since only even powers appear in the 
expansions of these combinations in powers of /9, they are analytic functions of 7, 

III. OTHER BOUNDARY CONDITIONS 

In the remaining portion of the paper we consider five different boundary condition at 
the ends of the polymer. In the case of two fixed ends with x = Xq = and f = fo = 0, the 
appropriate partition function is 

45(t)fixed,fixed = Zfl (0,0;0,0;t) = n-\bPf'^ [a'^ smh^ («t) - /3"2 sin^ {pi)y"\ (20) 

as follows from Eq. ffT9l) . For one free and one similarly fixed end, the partition function is 

/oo roQ 

dx / dvZfl (x,w;0,0;t) 
00 J —00 

= 2 [(1 + 27)a-2 sinh' (at ) - (1 - 2^) (3'^ sin^ (/3t ) + 4] "'^' . (21) 
In the case of two free ends, most relevant to experiments, 

/oo /"OO /"OO /"OO 

dx dv dxo dvoZy^l{x,v;xo,vo;t) 
00 J —00 J —00 J ~oo 

= 47r(6P)-i/2 [(1 + 27)2^-2 sinh" [at] - (1 - 27)^/3-2 sin^ (/3t )]"'/' .(22) 
For periodic boundary conditions 

^f6(^)periodic = / dxl dvl Z^^ {x,v; X, v;t) 



^'/,6 ^^)p«^"°'ii'= ~ / dx dv Zjf^ 

-00 J —00 J —00 



= - [cosh(at) -cos(/3t)] \ (23) 

In calculating the radial distribution function, Levi and Mecke p^] allow the polymer 
end points to fluctuate but with the end slopes fixed to the value zero. The corresponding 
partition function is given by 

/oo /»oo /»oo 

dx / dxo / Zfj (x, 0; Xq, 0; t) 
00 J —00 J —00 

= b-'/^P'/' [cosh^ (at) - cos^ (/3t )] "'/' (24) 

We have checked that Eq. ( p4|) is in complete agreement with the result |19| 



.(2) 



y(^) /,\ 00 

^/^bl'-Jlm n=l 



P{mi/ty + f{mx/tf + h 
Pin-n/ty + (s + f){mi/ty + 6 

6 



1/2 



(25) 



obtained by substituting x(t') = ^^gX„ cos(7mt'/t), < t' < t in the Hamiltonian and 

then performing Gaussian integrals over the coefficients Xn to obtain the partition function. 

For periodic boundaries the appropriate Fourier expansion is x{t') = ^^_oo^n^^''*"* ^*- 

This leads to 

zSf,,it),,r ^ fr P{2n7r/t)' + f{2n7r/tr + b 
45lWper l{P{2nn/tr + is + f)i2nn/ty + b' ^ ^ 

which is completely consistent with our result ( 123|) . 

For all of the boundary conditions considered above, Z\f^ has the asymptotic behavior 

-ln45(t)^«t = Eo(/,6K (27) 

^o(/, b) = a ( A^ '^' = 2-V2 (i) '^' [1 + i /(6P)-V2] V2 (28) 

both in the long-polymer regime t — )■ oo with / fixed and in the strong-stretching regime 
/ — )■ oo with t finite. In the latter regime Eo{f, h) in Eq. (!28|) further simplifies to E^i^f, h) ~ 






I (//P)^/^. As discussed in the appendix of Ref. [l7|, the dominant contribution to Zjl for 
large t comes from the "ground state" of the t-independent Fokker-Planck equation, which 
has eigenvalue EQ{f,b). 



IV. MOMENTS AND RADIAL DISTRIBUTION FUNCTION 

According to Eqs. (j4]) and ([5]), the mean value and variance of the difference ^ = L — Py = 
I Jq dt v(tY of the contour and longitudinal lengths are given by 

<^) = -^i^f41w4iwl (29) 



df 



d^ 



,(2) /.w(2) 



((e-(or) = ^in[4iw^}:4WJ (30) 

In Figs. 1 and 2, these two moments are plotted as functions of the contour length t = L 
and the longitudinal force parameter / = r/ksT for all five boundary conditions considered 
above. In both figures bx = by = 262P^^. As explained in footnote 2^, these are appropriate 
potential parameters for a long polymer confined in a channel with a square D x D cross 
section with D = ^P. 

That the lowermost and uppermost curves in Figs. 1 and 2 correspond to fixed-fixed and 
free-free boundary conditions at the ends of the polymer is reasonable, since these boundary 
conditions most and least restrict the polymer fluctuations, respectively. 



The leading asymptotic forms for (^) and ((^ — (0) ) fo^ large t and/or large /, toward 
which the curves in Fig. 1 and 2 tend, are the same for all five boundary conditions and 
given by 

iO ^ 2-'/'p-^" [h-'l' [I + \ fiKPyy^] -"^ + (6. ^ 6,)} t, (31) 

((e - (0)') ^ 2-9/2p~^/^ [h-'l' [1 + i f{KP)-'l\'"" + (6. ^ 6,)} t, (32) 

as follows from Eqs. (p7j) - (l30|) . In the strong-stretching regime / — )■ oo, Eqs. (13T|) and (132|) 
further simplify to 

(0 ^ \ ifPy/'t, ((e - (0)') ^ I /-=^/^P-^/^t, (33) 

independent of the potential parameters bx and by, as expected, and in agreement with the 



result of Marco and Siggia 23| for a strongly stretched, unconfined semiflexible polymer. 



From Eqs. ([T]) and f l3T|) we see that 

/* dt' v{t'f ^ {8bl/'P'/' + 4/P) "'/' + (6,. ^ 6,) (34) 



^2 



t 

for large t. From this relation one can identify the range of parameters b^, by, P, and / 
consistent with the domain of validity f ^ <^ 1 of our results. For b^ = by = 262P~^ and 
/ = 0, as considered in Figs. 1-3, Eq. fjM|) yields f^ = 0.176. 

We now turn to the probability distribution or radial distribution function 

m) = {s[^-l£dtv{tr^^ , (35) 

from which the above moments follow. Its Laplace transform is given by 



Pfis) = J^°^d^e-^^Pj{0 = (^exp(^-'-l^ dt', 



v{ty . (36) 



Comparing this expression with Eqs. (jlj) and (jS]), we see that 



41 (^) 41 (^) 



^n"") - ^(2) ,,, ^(2) ... ■ ^37) 



Substituting the partition functions (120|) -( !24|) in Eq. ( 137|) . we obtain simple analytic 
expressions for Pf{s) for each of the five boundary conditions considered above. The various 
moments of i^ may be generated from Pf{s) according to 

{C)={-d/dsrPf{s) , {{^-{Or) = i-d/dsr\nPf{s) , (38) 

s=0 s=0 



as follows from Eq. (136|) . and the radial distribution function -P/(0 is determined by the 
inverse Laplace transformation 

-| pc+ioo 

^/(0 = ^y dsPf{s)e^\ (39) 

Since -P/(s), as given by Eq. fl371) . only depends on s in the combination s + f, Eq. fj39|) 
implies 

PfiO = Aje-f^PoiO- (40) 

Thus, the radial distributions in the presence and absence of the stretching force only differ 
by a factor e~^^ and a corresponding normalization constant Af. From Eq. 036p and the 
normalization J^ dC, Pf{^) = 1, we see that Af = Po(/)~^; iii terms of the quantity Pf{s) 
defined in Eq. (|3j 



For the boundary condition that the polymer slope vanish at the end points and for / = 
and bx = by = b, Levi and Mecke [19'] obtained an analytical expression in the form of a 
power series for the radial distribution function P/(0 by substituting Eq. ( I25l) in Eq. (1371) 
and performing the integration over s in Eq. ( l39l) using Cauchy's residue theorem. For the 
periodic boundary condition Pf{C,) can be evaluated analytically, beginning with Eq. fl26l) . 
in the same way. 

For all five boundary conditions, for L = P and / = 0, and for the same potential 
parameters bx = by = 262P~^ as in Figs. 1 and 2, we have evaluated PjiO by performing 
the inverse Laplace transform ( l39l) numerically using Stehfest's method [24]. This leads to 
then distributions of R\\/L shown in Fig. 3. As the polymer is fairly short and not so 
tightly confined, the boundary conditions significantly affect the results. As expected, the 
narrowest, rightmost peak and the broadest, leftmost peak correspond to fixed-fixed and 
free-free boundary conditions, respectively. 

V. CONCLUDING REMARKS 

The properties of a semiflexible polymer confined along a line by a parabolic potential have 
already been studied in several other papers [M, [l6|420|, and we conclude by summarizing 



what this paper adds. As in Ref . p^ , we concentrate on the regime in which the polymer is 
stretched out along the line due to a strong confining potential and/or a strong longitudinal 
force applied to the ends. The main new result is an exact analytical expression for the 



partition function or path integral Z (r, v\ ro,Vo;t) , defined by Eqs. @ and (|5]), where the 
Hamiltonian (|2]) is quadratic in all three quantities df/dt, df/dt and r. Inclusion of the term 
quadratic in df/dt, with coefficient /, allows us to study the response of the polymer to a 
longitudinal stretching force and to calculate the distribution of the end-to-end distance or 
radial distribution function by replacing / by / + s in the partition function and performing 
the inverse Laplace transformation s — )■ ^ = L — i?||, as in Eqs. ( 137|1 and (139|) . 

From our result for the partition function Z (r, v; tq, vq] t) for arbitrary fixed endpoints 
and endslopes, we obtain compact analytical expressions (!20|) -( l24|) for the partition functions 
corresponding to five other boundary conditions of interest, including the case of freely 
fiuctuating ends most relevant to experiments. From these expressions it is a simple to 
calculate the radial distribution function by numerical inversion of the Laplace transform. 
In Fig. 1-3 we present numerical results, obtained this way, for the radial distribution 
function and its first two moments, for polymers which are short enough [L = P = 3D), so 
that the effects of the different boundary conditions at the ends are appreciable. 
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FIG. 1: Dependence of {^) = L — {R\\), where L is the contour length and Ru is the longitudinal 
length of the polymer, for fixed persistence length P, on L (upper figure) and on the longitudinal 
force parameter / = r/kBT (lower figure). As explained following Eq. (]30p . the potential parame- 



ters bx 



262P '^ used for the figures are appropriate for a polymer in a channel with a square 



D X D cross section with D = ^P. In the upper figure / = 0, and in the lower figure L = P. The 

curves are labeled according to the boundary condition at the ends of the polymer: free- free (rr), 

free-fixed (rx), Levi-Mecke (Im), periodic (pr), fixed- fixed (xx). 
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FIG. 2: Same as Fig. 1 except that the variance (A^^) = ((.^ — (,^)) ) is shown. 
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FIG. 3: Distribution of R\\/L for a polymer with contour length L = P with longitudinal force 



parameter / = 0. As in Figs. 1 and 2, the potential parameters are h^ 



262P ^, appropriate 



for a long polymer in a channel with a square D x D cross section with D = -^P. The curves are 
labeled according to the boundary condition at the ends of the polymer: free-free (rr), free-fixed 
(rx), Levi-Mecke (Im), periodic (pr), fixed- fixed (xx). 
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